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FREE-FLIGHT INVESTIGATION TO OBTAIN DRAG-AT-LIFT 
AND STABILITY DATA FOR A 60° DELTA-WING— BODY 
CONFIGURATION OVER A MACH NUMBER 
RANGE OF 1.5 TO 1.6 
By Clement J. Welsh 


SUMMARY 


Flight tests have been made of a rocket -propelled 60° delta-wing — 
body config\rrationj the wing airfoil section was NACA 0005-65 . Drag 
at lift and stability data were obtained for a Nfech number range 
of 1.28 to 1.6. 

The drag due to lift parameter increased from O. 5 I 'to O .58 with 
increasing Mach numbers, and only a small amount of leading-edge suction 
was obtained. Values of lift coefficient and pitching-moment coefficient 
varied linearly with angle of attack, and the aerodynamic-center position 
was nearly constant at approximately h6 percent of the mean aerodynamic 
chord. The damping-in-pitch derivatives varied from -1.6 to -1.2 with 
increasing Nfech numbers. Aeroelasticity effects calculated for the 
tested configuration increased with increasing Mach numbers and indi- 
cated a loss in lift of about I 5 percent and a forward movement of the 
aerodynamic center of approximately 4 percent of the mean aerodynamic 
chord at a Mach number of 1.6. 


INTRODUCTION 


As part of the National Advisory Committee for Aeronautics drag- 
due-to-lift program (see ref. 1, for example) a free-fli^t test has 
been made to furnish additional data for thin wings at high Reynolds 
numbers. A 6o° delta-wing configuration having an NACA 0005-65 airfoil 
section was flight tested and data obtained at Reynolds numbers around 

15 X 10^. 
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Also obtained and presented in this paper are zero-lift drag, lift, 
and stability data, including aeroelasticity corrections. The rocket- 
propelled-model configuration was tested at the Langley Pilotless Air- 
craft Research Station, Wallops Island, Va. 



SYMBOLS 

M 

Mach number 

q 

dynamic pressure 

R 

Reynolds number 

hf 

moment of inertia about the lateral axis 

s 

total wing area 

b 

wing span 

c 

local chord 

c 

mean aerodynamic chord 

L 

total lift 

D 

total drag 

M 

pitching moment 

*^L 

lift coefficient, L/qS 

^D 

drag coefficient, D/qS 

Cm 

pitching-moment coefficient (relative to the center of 
gravity) , M/qSc 

Cz 

section lift coefficient 

*^^pt 

optimum lift coefficient, corresponding to (L/^)max 

CDo 

zero-lift drag coefficient 

or K 

dCL2 

drag due to lift parameter 
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(Cmq + Cni(ij 


y 



0 

a 


lift-curve slope 
pitching -moment curve slope 

SCm 

damping-in-pitch derivatives^ — + 

Hi) 

lateral distance from fuselage center line 

pitching angle or local wing twist, in radians 
angle of attack, deg 



d = ^ da 

57.5 dt 

9 = de/dt 
Subscripts : 
r 

CP 

f 


s 


n-s 


rigid-wing conditions 

spanwise center of pressure 

flexible-wing conditions 

condition of full leading-edge suction 

condition of zero leading-edge suction 


MODEL AND TESTS 


A sketch of the model configuration is shown in figure 1, and a 
photograph in figure 2. 

The 60^ delta wing had an NACA OOO5-65 airfoil section and was con- 
structed of solid dural; the parabolic -profile fuselage was of wood- 
metal construction. The fuselage ordinates are listed in table I. 
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The model was equipped with a telemeter which transmitted measure- 
ments including angle of attack, free-stream total pressure, and normal, 
longitudinal, and transverse accelerations. 

Pulse rockets were installed in the rear portion of the fuselage 
to produce disturbances of the model in pitch about its lateral axis 
during flight in order that drag due to lift and stability data could 
be obtained. 


Preflight Measurements 

Experimental structural Influence coefficients were obtained for 
the 6o° delta wing corresponding to a loading along the 46-percent chord 
line. The coefficients which were used in the aeroelasticity corrections 
are listed in the following table as 9ij values corresponding to the 
local wing twist at station (i) due to a \mit load at station ( j) . For 
simplicity in listing these values, (i) and (j) are given integer val- 
ues 1, 2, 5^ 4, and 5 corresponding to actml spanwise stations (in 
inches) 4.40, 6.88, 9*10, 11.55^ and 14.48. The maxim\jm body radius 
is 5*25 inches and the actual semispan with the tip radius as shown in 
figure 1 is l6.75 inches. The values of (rad per lb) as listed 

are ^-10^] times the true values; 



During flight, the model was tracked with an ,NACA modified SCR-584 
radar unit to obtain position-time data and with a Doppler radar unit to 
obtain velocity-time data. Atmospheric conditions were determined from 
radiosonde equipment. The variations of wind direction and speed with 
altitude were obtained by tracking the ascending radiosonde with Rawin 
equipment. The variations of Reynolds number and dynamic pressure of 
the test with Mach number are shown in figure 5- 
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Acciaracy 

Experience has indicated that the errors in the absolute value of 
a measured telemeter quantity are within "tl percent of the range of the 
instr\imentj hence, at M = the errors in the normal and chordwise 

force coefficients have been calculated to be within ^ 0.005 and -0.001, 
respectively. Errors in Mach niomber and Ax are believed to be within 
iO.Ol and ±0.10°, respectively, throiighout the test range. 


Analysis 

The disturbances in pitch resulting from the pulse rockets resulted 
in a maximum angle of attack of about t6°. As there were no appreciable 
transverse accelerations during flight, the short-period oscillations 
resulting from these disturbances have been analyzed assuming two 
degrees of freedom by the methods of reference 2. In addition, the 
instantaneous pitching-moment coefficients have been calculated from 
the two normal accelerometers as indicated in reference 3* 

Expres sions for and (l/d)j^j„ = 

( i ) = -i / 1 are obtained from the general expression for drag 

\D/max ^^Dq 

coefficient, Cj) = + KC^^. The upper range of the test data 

was in the range of hence the parabolic variation of Cp 

against was well defined. 

To indicate the aeroelasticity effects on the lift and static 

longitudinal stability resulting from the flexibility of the 3 -percent- 
thick solid d\oral wing, a brief aeroelasticity analysis was made. The 
method used is indicated in the appendix and makes use of wing struc- 
t\iral influence coefficients combined with a spanwise loading for the 
rigid wing. As previously stated, the influence coefficients were 
obtained experimentally, while the spanwise loading used was that for 
a 68° delta wing at M = 1.6 (experimental curve from unpublished 
data) . The difference in loading between the 6o° and the 68° delta 
was assumed negligible. 



DISCUSSION 


The data in general are presented for a wing-body configuration 
having a wing with flexibility characteristics that could be represent- 
ative of that of a typical missile or aircraft in this speed range. 
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In specific cases, ciirves are shown which have been corrected for aero- 
elasticity to indicate the range of losses in lift and static longitudi- 
nal stability that could exist on similar configurations relative to a 
rigid-wing condition. A plot of the flexibility parameter, effective- 
lift ratio, for the configuration is shown in figiire 4(a) . 

To confirm that the assijmptions made in applying the method for 
aeroelasticity corrections indicated in the appendix to the delta-wing 
configuration of the present test were Justified, values of the 

present test corrected to rigid-wing conditions and two similar configu- 
rations (ref. 4) are shown in figure 4(b) . The two configurations frcan 
reference 4 furnish a good comparison being delta-wing — body configura- 
tions having wings with identical sections to that of the present test 
and mounted on similar bodies. As the wings were constructed of solid 
steel, they are assumed to be representative of nearly rigid wings. 


Drag 


In figiire 5 typical plots of Cq against and Cp against 

are shown, including test points, to give an indication of the scatter 
of the actual data points. 

In figures 6(a) and 6(b) , zero-lift drag and drag due to lift are 

d-C 

shown against Mach number. The values of — increased from 0.3I 

dCr^ 


at M = 1.28 to 0.58 at M = 1.6. Shown for comparison is a calculated 


cirrve of 



(no leading-edge suction for the flexible wing con- 


n-s 


figuration) . The 


^dC 


D 


dCt 


values were calculated from the integrated 


in-s 


effect over the wing resulting from the assumption that the resultant sec 
tion force was normal to its local chord and in consideration of the span 
wise varying a resulting from the flexibility of the wing. Also shown 

is a / — ciirve, representative of a condition of full leading-edge 




suction. The 


'dCr 


dCD^ 


curve was obtained by subtracting A 


dCy 


dCL2 


(the increment of the calculated leading-edge contribution to the 


values 


dCy 


dCi 
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of an isolated 6o° delta wing when full leading-edge suction exists, 

values . 

n-s 

As would be expected and as is apparent from these curves^ the 

dCn 

experimental curve of — ^ for the configuration includes a combined 

dCL^ 

effect of leading-edge suction and aeroelasticity^ and indicated that 
only a small amount of leading-edge suction was obtained. The curves 
also Indicate the major effect that leading-edge suction can have on 
the drag at lift of a 6o^ delta wing over this speed range. 

To give an indication of the performance characteristics of the 
configuration, (V^)max ^^opt been plotted against Mach 

number in figures 7 ( 0 -) Q'^d T(^) • (V^)max decreased from 6.6 to 5*8 
with increasing Mach number; remained nearly constant over the 

test range at approximately 0.24. 


ref. 5) from 


^dC 

dCr 


D 


Stability 


Both and values varied linearly with a throughout the 

test range and their corresponding slopes are plotted in figures 8(a) 
and 8(b) against Mach niomber; both C-^ and decreased with 


increasing Mach number. Also shown are 


Ci^ and 


Cm^ curves corrected 


for aeroelasticity. 


As is obvious, the aeroelasticity effects increased with increasing 
Mach numbers, and at M = 1.6 the loss in lift was about I5 percent 
and the loss in Crn^ was approximately 25 percent ^the loss in 

relative to a pitching axis at 14.6 percent MAC^. 

Calculated values are shown as determined by reference 5 fo^ 

an isolated 6o^ wing. Also shown are calculated wing -body Cj^ values 

as determined by reference 6. The agreement between calculated values 
and the corrected experimental values is very good. 

The aerodynsimic- center position is shown in figure 9{^) against 
Mach number and was. approximately constant over the test range at about 
46 percent MAC. With aeroelasticity corrections as shown, the aero- 
dynamic center moved rearward; at M = I.60 this change in the 
aerodynamic -center postion was about 4 percent MAC. 
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Shown in figure 9(^) the damping-in-pitch derivatives 4 - 

against Mach number. The sum of the derivatives varied from -1.6 to 
-1.2 for corresponding Mach numbers of M = 1.28 and M = 1.60. Cal- 
culated damping-in-pitch derivatives determined by reference 7 also 
shown; the excellent agreement with the measured values would appear 
fortuitous. Damping-in-pitch derivatives from tunnel measurements, 
reference 8, are also in agreement with the present test. 


CONCLUSIONS 


The results of the present investigation of a 6o° delta-wing — body 
configuration with an NACA 0005-65 airfoil section indicate the following 
conclusions for a Mach number range between M = 1.28 to M = 1.60: 

1. The drag-due-to-lift parameter increased from O.5I to O.58 with 
increasing Mach number, and only a small amount of leading-edge suction 
was obtained. 

2. Cl and C^^^ values varied linearly with a and the aerodynamic- 
center position was nearly constant at approximately 46 percent MAC. 

5- The damping-in-pitch derivatives varied from -1.6 to -1.2 with 
increasing Mach number. 

4. Values of (V^)max (decreased from 6.6 to 5*8; however, 

C^opt remained nearly constant at 0.24 over the Mach number range. 

5 . Aeroelasticity effects calculated for the tested configuration 
increased with increasing Mach numbers and indicated a loss in lift of 
about 15 percent and a forward movement of the aerodynamic center of 
approximately 4 percent MAC at M = 1.6. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., July 8, 1955- 
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APPENDIX 


AEROELASTICITY CORRECTIONS 


Aeroelasticity corrections to the wing-body lift and static longi- 
tudinal stability data have been made as was similarly made for exposed 
swept wings in reference 8, The derivation for flexibility effects for 
the wing-body configuration is presented as it is slightly different 
than when only the exposed wing panel is considered and in turn results 
in immediate wing-body flexibility effects. 

The basic assumptions made in making the present corrections are: 

(1) the spanwise loading of the rigid wing mounted on the body is the 
same as for the isolated wing, that is, the body effects are negligible, 

(2) that the loaded flexible wing has a straight chord line at any span- 
wise station, (3) that in determining the structural influence coeffi- 
cients, applying the loading along the approximate center of pressure 
chordline of the wing will result in coefficients representative of the 
loading in flight, and (h) aerodynamic induction effects are considered 
negligible after the initial rigid-wing loading. 

The method requires both the rigid-wing spanwise loading curve and 
structural influence coefficients, and is essentially a form of strip 
theory in that the wing and its corresponding loading curve is divided 
into a series of strips (may be of varying widths) similar to the sketches 
below, where the (O) strip corresponds to the portion of the semispan 
within the fuselage. 



Wing Panel Spanwise Wing Loading 
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From the previous assumptions it is apparent that the lift of each strip 
of the flexible wing is directly proportional to (a^ + Ax) where xp is 
the angle of attack at the center line of the configuration or a for 
rigid-wing conditions, and Zki due to flexibility is the incremental a 
at the center of pressure of the strip. As the equilibrium Ax for any 
strip is dependent on the twisting of the other strips, the resulting 
lift on a strip will he given by the solution of a series of equations 
equal in number to the N number of strips into which the wing is 
divided. 

The total lift over the semispan of a rigid wing is given as 
L = Cj;^(Xj.qS/2, where a is in radians and the subscript r refers to 

center line or rigid conditions. Hence, the lift associated with any 
X-strip of the flexible wing would be given by 

Lx = C]^^q(s/2)Kx(ar + ^x) 

where (from the loading curve) 


K 


X " 



C^c 

CTc 

^L^av 


dr\ 


By dividing the expression for Lx hy an expression for 

obtained 





o-r Ci^q(s/2)Kx 


1 


Ax Y 

— ^ is 

Xp 


and to simplify for computational purposes 

^ % 1 _ 

Oj. Oj. QK'x 


( 1 ) 


where 


K'x= KxS/2 

the number of equations resulting being equal to the N n\miber of strips 
chosen. 
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The incremental a due to flexibility can also be written in terms 
of the structural influence coefficients and for the Ax associated with 
the typical X-strip is, 


AXy L-| Lq Ly Lkt 

^ ®X1 + — 9X2 + . . . — exx + • . . - 0XN (2) 

where the general influence coefficient is defined as the twist 

per unit load at the center of pressure of the strip i due to a load 
at the center of pressure of strip j. Again the number of equations 
resulting is equal to the N number of strips chosen. 

Combining equations (l) and (2) and their corresponding expressions, 
yields the following set of equations: 


corresponding to the strip within the fuselage as indicated in the pre- 
vious sketches, and then the N-1 simultaneous equations: 



ar 




®12 + 



IX + 



®1N = 


-1 




9vi + ®YP 

AJ . (Xr ^ 


% 

(Xp 



e 

Xp 


XN - 


-1 


Li 


+ 




% 

Xp 



-1 


Lq L-j 

The above equations are solved for — , 

ax» ar 

in terms of specific values of Q (the values are chosen corresponding 
to the test range) . With the solution of the previous equations the 
following expression can be written. 


L2 Lx % 

, • • • “■ ", • • • 

Xp Xp Xp 





_x 


W2 


% 

Oj. 
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which is the ratio of the lift of the flexible-wing configuration to 


of the above expression for the assumed .specific values of Q previously 
referred to, for any given q value and in turn its corresponding Mach 
number . 

The spanwise center of pressure for the flexible wing is given by 


where A is the approximate sweep of the center -of -pres sure chord line. 


As previously mentioned, the method requires a linear chordwise 
variation of the flexible-wing chord plane at any spanwise station} 
hence, it seemed questionable that it could be used for delta wings in 
general. However, in finding the experimental influence coefficients, 
it was found that with a spanwise loading along the 46 percent chord 
line (representative of 60 ° delta wings in this speed range) resulted 
in a linear variation for this wing within the accuracy of the measure- 
ments. The final justification of assuming a linear chordwise variation 
is in the agreement of the corrected lift values with those measured 
under near rigid-wing conditions. 

For the present test, the wing was divided into six strips which, 
from previous tests, seemed to be an adequate number. Though inertia 
effects can easily be included in the present method, a rough check 
indicated the inertia effects of the present wing were negligible} hence, 
no corrections for these effects were included. 




hence 



L 


And finally is given by 
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TABLE I 


FUSELAGE ORDINATES 


Distance from nose 
of fuselage, in. 

Fuselage radius, in. 

0 

0 

.590 

.097 

.585 

.145 

.975 

.259 

1.950 

.469 

5.900 

.902 

5.850 

1.298 

7.800 

1.658 

11.700 

2.267 

15.600 

2.730 

19.500 

5.047 

25.400 

5.218 

27.500 

5.248 

51.200 

5.221 

55.100 

5.161 

59.000 

5.069 

42.900 

2.945 

46.800 

2.785 

50.700 

2.594 

54 . 600 

2.571 

58.500 

2.115 

62.400 

1.826 

65.000 

1.615 
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Figure 1.- Sketch of the tested configuration, all linear dimensions 

in inches. 
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Figure 2.- Photograph of the tested configuration. 
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Figure 5-- Variations of Reynolds number based on mean aerodynamic chord 

and q with Mach number. 


NACA RM L55G14 CONFIDENTIAL 


18 


CONFIDENTIAL 


NACA RM L55G1^ 



Figure A.- Lift characteristics of the tested configuration and con^jari- 

sons with near rigid wings. 
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Figtire 5*- Typical Cjj curves showing actual test points. (First cycle 

of M = 1.5 oscillation.) 
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Figure 6.- Variations of zero-lift drag and drag at lift with Mach niomber. 
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Figure J.- Variations of maximum lift-drag ratio and optim\am lift 

coefficient with Mach number. 
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Figure 8.- Variations of the slopes of the lift and moment curves with 
Mach number including aeroelasticity corrections. 
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Figure 9-~ Variations of the aerodynamic center and damping in pitch 

derivatives with Mach mmiber. 
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